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We have used a superconducting single-electron transistor as a DC-electrometer that is strongly 
coupled to the metal island of another transistor. With this set-up, it is possible to directly measure 
the charge distribution on this island. The strong capacitive coupling was achieved by a multilayer 
fabrication technique that allowed us to make the coupling capacitance bigger than the junction ca- 
pacitances. Simulations of this system were done using orthodox theory of single-electron tunnelling 
and showed excellent agreement with the measurements. 



A single-electron transistor can make extremely sensitive charge measurements with a resolution of about 
1CP 5 e/VHzJjl3 In principle, this sensitivity is sufficient to measure the charging and discharging of a small con- 
ductor as current flows through that conductor. A single-electron transistor would be able to register the tunnelling 
of individual electrons tunnel onto the conductor from one lead and off the conductor to another lead for currents up 
to a few picoamperes. In practice, the high output impedance of a single-electron transistor makes this measurement 
very difficult. If the charge on the conductor as a function of time could be measured, then it would be possible to 
determine the occupation probabilities of the various charge states. In charge state (0) the conductor is neutrally 
charged, in charge state (1) the conductor has an excess charge of one electron, and so forth. The charge state 
occupation probabilities play a central role in the theory of single-electron tunnellingEJLl but are typically not directly 
accessible experimentally. Here we report on an experiment where a superconducting single-electron transistor (SSET) 
was used to directly measure the charge state occupation probabilities on the island of another nearby SSET. 

An essential feature of this experiment was that the measurement SSET (the electrometer) was strongly coupled 
to the island where the current was flowing through. Strong coupling means that the coupling capacitance was 
comparable to the total capacitance of the nearby island. When one electron was added to this island, about e/4 of 
charge was induced on the coupling capacitor. The electrometer was biased at a small voltage and it's tuning gate 
was used to scan the charge on the neighboring island. These gate traces directly reflect the charge distribution on 
this island. The measurements are consistent with orthodox theory and they show that a strongly coupled SSET can 
be used to directly measure the charge distribution. 

I. EXPERIMENT 

A schematic of the device is shown in Fig. |l|a. The device was fabricated in three layers. The junctions were 
fabricated using standard shadow evaporation of aluminum. SSET1 has a planar gate capacitor C g \ while the gate 
capacitor G g 2 is defined as a parallel plate capacitor. Details of the fabrication of a similar device were described 
elsewhere.Q Figure |l|b shows a SEM picture of the device. The two square islands of the SSETs are coupled via an 
underlying dumbbell shaped coupling capacitor. The total effective capacitance between the two islands is called 
C m . Both SSETs were biased asymmetrically, connected to a voltage source at one side and grounded at the other 
side. The device was measured in a dilution refrigerator with a base temperature of 10 mK. The leads were equipped 
with 7r-filters at room temperature and standard copper-powder filters at base temperature. The effective electron 
temperature, 25 mK, was measured in the normal state by fitting experimentally obtained Coulomb peaks. All further 
measurements were done in the superconducting state, the superconducting gap being A = 200 /ieV. 

Throughout the measurements, the voltage bias of the electrometer was kept constant at Hi = 805 fj,V, just above 
4A/e. The current through the electrometer I\ was measured as a function of the gate voltage V g i and the bias 
voltage Vb2 of the SSET2. The gate voltage V g 2 was kept constant. Figure |]a shows typical Coulomb oscillations 
of the current through the electrometer. The gate of the electrometer was swept while SSET2 was biased at 800 
^V, the current I2 being negligibly small. Figure IJd shows the same Coulomb trace when Vb2 — 890 ^V, above the 
quasiparticle threshold of SSET2. Surprisingly, the Coulomb peak is split into two peaks, while at even higher bias 
{Vb2 = 1090 /iV) it is split into three. As we will explain below, each extra peak can be attributed to the presence of 
an extra electron on the second island. 
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When Vb2 = 800 /iV, the current through SSET2 is still negligibly small and the occupation probability of charge 
state (0) on island 2 is nearly 1. When Vb 2 is higher than the threshold voltage, a quasiparticle current I2 will start 
to flow and the charge 712 on the island 2 will switch between (0) and (1). The presence of an extra electron on island 
2 will induce a fraction of an electron on coupling capacitor C m . By writing down the total charge on both islands 
as a function of the capacitances and the island potentials, one can show that this fraction is C m /C-£2, where Cy/i is 
the sum of all capacitors connected directly to island 2. Table | gives the capacitance and resistance values of all the 
circuit elements. With this table, we can calculate that the charge induced on capacitor C m is 0.27 e. This results 
in the extra Coulomb peak (labelled (1)) shifted —e/C g i * 0.27 = —1.04 mV with respect to the peak labelled (0) 
in Fij* ||b. When Vfo is increased even more, charge state (2) is also populated on island 2 and three peaks appear 
(Fig. gc) . Figure [| shows the Coulomb traces of the electrometer current I\ in grayscale versus the bias Vfo ■ One can 
clearly see that consecutive charge states become populated with increasing V^- The average bias voltage difference 
between successive charge states on island 2 is 2Ec2/e, where Eq2 is the charging energy of island 2. The charge state 
(4) becomes populated at Vb2 = 1180 /iV and induces 1.08 e on island 1. The corresponding peak in the Coulomb 
trace of the electrometer overlaps with the next set of Coulomb peaks, limiting the number of observable charge states 
to four with these circuit parameters. 

By measuring Coulomb oscillations, the electrometer can be used to directly resolve the average population of charge 
states on a nearby island, even though the charge on this island changes on a nanosecond timescale. At Coulomb 
peak (j), current can only flow through the electrometer when the island 2 is in charge state (j). The fraction of the 
time that island 2 spends in charge state (j) is equal to the relative peak height defined by, 

where 7 p j is the peak height of the Coulomb peak corresponding to the charge state (j) on island 2. Simulations 
confirm that the height of the individual peaks reflects the exact population of the corresponding charge state. 

The relative peak heights of the Coulomb traces in Fig. ^b,c are calculated and shown in column I of Table |l[ They 
closely match the occupation of the various charge states as calculated in the simulations for identical bias conditions, 
shown in column II. Also, for bias conditions other than Vb2 = 890 fiV and 1020 fiV, the simulated relative peak 
heights closely match the experimental ones. This shows that a strongly coupled SSET can be used to quantitatively 
measure the charge distribution on a nearby object. 

II. SIMULATIONS 

The current through both SSETs was calculated using a master equation analysis. By solving the master equation 
we can calculate the occupation probabilities of the various charge states. The master equation for this two island 
system is, 

= ^ {Pkl^kl^ij — Pij^ij-t-kl) ) (2) 
kl^ij 

ij 

where Py is the probability that the system has i excess electrons on island 1 and j excess electrons on island 2. 
r denotes the transition rate between different charge states. In the stationary state, Py does not change and the 
left hand side of equation |^ is zero. The first term on the right describes the population of charge state ij from 
charge state kl while the second term describes the depopulation of charge state ij to charge state kl. We neglect 
co-tunnelling processes and T only is non-zero when either i — k±lorj = l±l. Furthermore we only take into 
account a maximum of five charge states per island. 

The superconducting tunnel rates T were then determined with Fecmi's Golden Rule using the superconducting 
density of states and the free energy difference AP of a tunnelling event Eu AP is the sum of the change in electrostatic 
energy plus the work done by the voltage sources. The total electrostatic energy of the system can be written as: 

P(ni,n 2 ) = E C i(n i + nx) 2 + E C 2{riQ2 + n 2 ) 2 + E m {n m + ni)(n 02 + n 2 ), (4) 

e 2 C 2 s 

E ci = 0(r n ^2"T = 95 ^eV, (5) 
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e 2 C 

Ec2 = W^c^cT) = 53 ^ (6) 

E m = f° m 2 = 51 M eV, (7) 

where tiqi and riQ2 are the normalized charges induced on the islands by the voltage sources, n\ and rii are the excess 
number of electrons on the islands and Cjji and Cs2 are the sum of all capacitors directly connected to the respective 
islands. The electrostatic energy has three contributions, the charging energies of the separate SSETs (equations ^|and 
6) and the coupling energy E m , which describes the electrostatic interaction between both SSETs. When equations 
2 and [| are solved for P^-, the total current I\ can be calculated with P^ and the tunnel rates. 

In Fig. H the current through the electrometer has been calculated for the same bias conditions as Fig. ^. One can 
clearly see the extra Coulomb peaks appear when the bias voltage Vb2 is increased. The absolute peak height of the 
experiments is about 60% of the peak height in the simulations. This can be accounted for by the rounding of the 
superconducting gap. Instead of the discontinuous jump in quasiparticle current through a superconducting junction 
at 2A/e, in real experiments, the current increases with a non-zero slope. In these experiments, the differential 
resistance in this regime is about 5% of the high bias junction resistance. For a bias voltage of Vbi = 805 /aV, only 
5 fiV above 4A/e, this has two consequences. First, the Coulomb peaks have a more triangular form as can be seen 
in Fig. ^], second, the Coulomb peak height is smaller than in the simulations where the rounding has not been taken 
into account. Simulations where the rounding of the gap was taken into account with a simple model showed that 
the rounding of the gap does not change the relative height of the peaks, it merely decreases the overall current . 

For the simulated Coulomb traces of Fig. [|, the relative peak heights as specified by Equation [l] are given in Table 0. 
The experiments closely match the simulated values. The relative peak heights in simulations arc slightly different 
though from the occupation of the charge states on island 2 when the electrometer is switched "off" (Vb2 = 800 /^V). 
Column three of Table |l| shows the undisturbed occupancies of the three charge states as determined from the 
population matrix Py. As can be seen from Table ||, the bias of the electrometer has a small back action on the 
occupation of charge states on island 2. For the bias range of Fig. |[ it can be shown that the back action of this 
electrometer changes the occupancies of the various charge states by a maximum of 5%. 

Figure ||is the simulated equivalent of Fig. [|. The extra peaks appear in Fig. || at exactly the same bias conditions as 
in Fig. |^, demonstrating the close agreement between experiments and simulations. Another feature that clearly shows 
up in the simulations as well as the measurements is the existence of a current plateau in between the neighboring 
Coulomb peaks. Under these bias conditions, the electron-tunnelling^tbrough both SSETs is correlated. This effect 
has been described before in coupled ID arrays of tunnel junctions. 0113 The details of this effect will be discussed 
below. 



III. DISCUSSION 



The ability to determine the position and the height of the extra Coulomb peaks gives constraints on the bias 
conditions. In general, the width of the peaks has to be smaller than the separation between adjacent peaks. In 
the superconducting state the width of Coulomb peaks is almost independent of temperature for ksT < 0.5A and 
depends linearly on the applied bias. This constraint can be rewritten as: 

eVu - 4 A < E m . (8) 

This simply states that the energy associated with the voltage bias has to be smaller than the coupling energy. Because 
of the quasiparticle threshold at 4A/e, this constraints the bias voltage to 800 fiV < V&i < 851 /iV for this sample. 
The quasiparticle rate is almost independent of the bias in this bias window and simulations indicate that the back 
action of the electrometer is also constant. If we take into account the rounding of the gap and the experimental 
current noise, Vb2 = 805 fiY is about the optimal bias voltage, combining an acceptable signal to noise ratio with a 
reasonably small width of the Coulomb peaks. With the current sample parameters, we are limited to the observation 
of a maximum of four charge states on the neighboring island. We estimate that it is feasible to observe at least 
seven different charge states, when the coupling capacitance is lowered to 190 aF, while keeping the other sample 
parameters constant. 

Both Fig. |^ and ^ clearly show the existence of a current plateau in between the accompanying Coulomb peaks. 
In order to be able to measure the relative peak heights, this plateau current should not exceed the Coulomb peak 
current and therefore its mechanism should be understood. The mechanism can be most easily explained when the 
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number of occupied electron states on island 2 is limited to two and under the assumption that the tunnel rates in 
SSET1 are much larger than those in SSET2. Figure ^a schematically displays the quasiparticle thresholds for SSET1. 
The position of the dots denotes the effective background charge when the charge state of island 2 is (0) (right dot) 
and (1) (left dot). The position of the dots relative to each other is fixed. The bias voltage Vbi and hence the dots lie 
just above 4A/e. With the gate voltage V g \ the position of both dots can be shifted along the Qoi axis. If the gate 
voltage positions one of the two dots above both quasiparticle thresholds a and (3 this leads to current in the form of 
a Coulomb peak. If the dots are positioned as depicted in Fig. ^a there is an additional mechanism that will carry 
current. 

The charge states with the lowest energy are now (10) and (01). If a current is forced to flow through SSET 2 by 
biasing it above its quasiparticle threshold, the following current cycle is most probable: If we start with charge state 
(00) it is only favorable for electrons to tunnel onto islands 1 or 2 via the top junctions. Because we assume that the 
tunnel rates in SSET1 are much larger than those in SSET2, an electron will most probably tunnel through junction 
jl first, as shown in Fig. ^p. Now the system is in the charge state (10) which is stable for electron tunnelling in 
SSET1. After some time the bias voltage Vbi forces an electron on island 2 and the system is in state (11). As can 
be seen in Fig. ^c this state decays to (01) through junction j2, again assuming the tunnel rates are much higher in 
SSET1. This charge state is also stable for electron tunnelling in SSET1. The cycle is completed when the electron 
is forced off island 2 and the system is back in charge state (00) . The cycle of one electron tunnelling through SSET2 
has transported another electron through SSET1 making I\ = Ii. This cycle is possible for all gate voltages where 
both the position of ni = (1) lies below quasiparticle threshold a and position of n 2 = (0) lies below quasiparticle 
threshold (3. This gives rise to a current plateau exactly in between the Coulomb peaks attributed to the both charge 
states. 

In this sample, the resistances of SSET1 and SSET2 and hence the tunnel rates differ by only a factor of two. This 
means that cycles can be missed, for example if state (00) decays to (01), the system is forced to (00). An electron is 
transported through SSET2, without giving rise to current in SSET1. The general equation for the relation between 
Ji and I2 can be deduced by analytically solving the master equation under the assumption that only the four charge 
states (00), (10), (01) and (11) need to be considered. If we assume that the tunnel rates through junction jl and j2 
are equal and called Ti, as well as those through j3 and j4 are equal and called T 2 this yields: 

By deriving the expressions for the Coulomb peak current, it can be shown that the peak currents are always sufficiently 
larger than the plateau current, making it possible to adequately determine the relative peak heights. When the 
number of occupied charge states on island 2 is larger than two, the mechanism leading to the current plateaus is 
similar, but different combinations of charge states might be stable and equation ^| will be modified. Again though, 
the plateau current is always smaller than the Coulomb peaks adjacent to the particular current plateau. 

We also studied the performance of the electrometer in the normal state. In the normal state, however, the 
Coulomb peaks are very sensitive to thermal fluctuations. The thermal broadening of the Coulomb peaks at 70 mK 
was enough to merge the adjacent Coulomb peaks, making an accurate determination of the relative peak heights 
impossible. Additionally, due to a mechanism similar to the one leading to plateau current in the superconducting 
state, the adjacent Coulomb peaks merged at 30 mK in the normal state, making the normal state operation of this 
electrometer impractical. 



IV. CONCLUSIONS 



We have used a SSET to measure the charge distribution on a neighboring island. Both the islands were strongly 
coupled by a multi-layer technology. The presence of an extra electron on a neighboring island split the Coulomb 
peaks of the SSET. The relative height of these peaks directly translates to the occupation of the associated charge 
state. In between the neighboring Coulomb peaks the current is carried by correlated tunnelling of electrons through 
both SSETs. 
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FIG. 1. (a) Schematic of the circuit. The two junctions on the left form the electrometer (SSET1), whose island is coupled 
capacitively to the island of a nearby SSET (SSET2). (b) Scanning electron microscope picture of the completed device. The 
light gray layer is fabricated in gold, the aluminum layer shows up as dark grey. 
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FIG. 2. Experimental Coulomb traces of the electrometer for different values of Vbi while Vbi = 5/iV, V g i = 10 yN and T=25 
mK. The extra peaks Coulomb peaks in (b) and (c) correspond to the presence of extra electrons on island 2. 
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FIG. 3. Electrometer current versus V g i and Vbi- White indicates no current, black indicates a maximum current of 25 pA. 
The arrows indicate the values of Vbi where the traces of Fig. ^a-c have been extracted. At Vbi = 1180 fiV the charge state (4) 
becomes populated, but the corresponding peak overlaps with the neighboring set of Coulomb peaks. 
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FIG. 4. Simulations of Coulomb traces of the electrometer for different values of Vb2- Hi = 805 /iV, V g i = 10 /iV and T=25 
mK. 
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FIG. 5. Simulation of the electrometer current versus V g i and Vb2- White indicates no current, black indicates a current of 
45 pA. The arrows indicate the values of Vm where the traces of Fig. 0a-c have been extracted. 
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FIG. 6. (a) Schematic of the quasiparticle thresholds of SSET1 above Hi = 4A/e, shown as thick lines. The positions 
of the dots denote the effective background charge induced by the absence or presence of an extra electron on island 2. (b) 
Visualization of the possible tunnelling events on island 1. When the system is in charge state (00) it will decay to (10) by an 
electron tunnelling through junction jl. Electron tunnelling through junction j2 is energetically unfavorable, just like electrons 
tunnelling upward, (c) In a similar way charge state (11) decays to (01). 
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TABLE I. Capacitance and resistance values for the circuit parameters as calculated from the stability diagrams and cur- 
rent-voltage characteristics of both SSETs. 
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TABLE II. The population of first three charge states on island 2, as calculated from the peak heights in the experiments 
(Fig. |^) and the simulations (Fig. ^) . The undisturbed population is determined by a calculation of the population matrix Pij 
of equation m when the electrometer is switched off (14i = 800/iV). 
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